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ABSTRACT 

NGC 5194 (M51a) is a grand-design spiral galaxy and undergoing interactions with 
its companion. Here we focus on investigating main properties of its star-formation 
history (SFH) by constructing a simple evolution model, which assumes that the disc 
builds up gradually by cold gas infall and the gas infall rate can be parameterizedly 
described by a Gaussian form. By comparing model predictions with the observed 
data, we discuss the probable range for free parameter in the model and then know 
more about the main properties of the evolution and SFH of M51a. We find that the 
model predictions are very sensitive to the free parameter and the model adopting a 
constant infall-peak time t p = 7.0Gyr can reproduce most of the observed constraints 
of M51a. Although our model does not assume the gas infall time-scale of the inner 
disc is shorter than that of the outer disc, our model predictions still show that the 
disc of M51a forms inside-out. We find that the mean stellar age of M51a is younger 
than that of the Milky Way, but older than that of the gas-rich disc galaxy UGC 8802. 
In this paper, we also introduce a ’toy’ model to allow an additional cold gas infall 
occurred recently to imitate the influence of the interaction between M51a and its 
companion. Our results show that the current molecular gas surface density, the SFR 
and the UV-band surface brightness are important quantities to trace the effects of 
recent interaction on galactic SF process. 

Key words: galaxies: abundance — galaxies: photometry — galaxies: evolution 
galaxies: individual: M51a — galaxies: spiral 


1 INTRODUCTION 

M51a is one of the nearest galaxies with grand-design spi¬ 
ral arms. The outstanding point of this galaxy is that 
it is known to be interacting with its companion galaxy 
NGC 5195 (the interacting system NGC 5194 + NGC 5195 
is often called as M51). Coupled with its large angular size, 
nearly face-on inclination, high surface brightness and kine¬ 
matic complexity, many studies about M51a have been car¬ 
ried out using various wavelengths and methods. 

With respect to dynamic study, M51a is a very good 
target to explore the nature and origin of spiral structure. 
Two opposite scenarios dominate the discussion in the liter¬ 
ature. One theory considers that the spiral arms are a long 
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lasting pattern that slowly evolv es and rotates with a single 
angular speed dLin fe Shulll96j f. while the other assumes 
the arms to be transient disturbances generated, e.g. , by 
the t idal interaction with a companion ( Toomre fe Toomrel 
1 1972h . In the seminal kinematic study of M51a, Tullvl ( 1974 1 
identified the spiral pattern in the outer disc as transient 
feature, while the inner arms are th ought to be in a steady 
state. Indeed, IColombo et al.l d2014l 'l showed the kinematic 
evidence of an m = 3 wave in the inner disc of M51a, 
which may suggest the density-wave nature of the main spi¬ 
ral structure of M51a. On the other hand, many kinematic 
and hydrodynamical models showed evidence that, instead 
of quasi-steady density wave, the internal structure of M51a 
drives from the complicated and dynamical interaction with 
its companion dD obbs et alJ l201(J : Salo fe Laurikainenll2000l : 
iTheis fe Spinnekeidl2003l l. 


In addition to its dynamics, M51a is also a test bed 
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to understand the effects of galaxy interacting on its star 
formation (SF) process by investigating its stellar popula¬ 
tions. Early observations have already shown that the blue 
and yellow-red stellar popul ations within M51 are spatially 
separate structures ( Zwick vl 119551 ). Later investigations of 
star clusters in M51a revealed that the star cluster forma¬ 
tion rate increased significantly during the period of 100- 
250 Myr ago, which is consistent with the e£och of the dy¬ 
namic encounters of thes e two galaxies (I L ee et ah l 20051: 
iHwang fe Leel 120081 . l2010h . iKaleida fe Scowen ( 2010l) also 
found that there was an enhancement in the number of 
young stellar associations in the northern arm closest to 
the companion, most likely triggered by the interaction of 
M51a and NGC 5195. Through the spectral energy distri¬ 
bution modeling of the m ulti-bands images of the M51 sys¬ 
tem, iMentuch Cooper et alj (120121') revealed a burst of star 
formation occurred in both galaxies roughly 340-500 Myr 
ago, which is in agreement with the result s of the colour- 
magn itude diagr ams of individu al stars dTikhonov et al.l 
l2009h . Moreover, iLee et ahl (l201ll ) found evidence that the 
tidal interaction between M51a and its companion appeared 
to enhance the SF process at the tidal bridge connecting the 
two galaxies. 

On the other hand, although the parameterized mod¬ 
els have already been applied to several individual galaxies 
and proven to be f ruitful tools to explore galactic f orma¬ 
tion and evolution dTinslevI 1980^_Changet_ak[ 199fjj_20121 


Boissier&Prantzos 20001: IChiappini et al.ll200ll ; lYin et aid 


20091 : Kang et alj 201^ ). there is still lack of similar inves¬ 


tigation on M51a partially due to the complexity caused 
by the on-going interacting between M51a and its compan¬ 
ion. In this paper, we construct a simple evolution model 
for M51a to build a bridge between its SFH and its ob¬ 
served properties, especially the radial distributions of cold 
gas surface density, metallicity, and the radial profiles of sur¬ 
face brightness in multi-bands. We present a parameterized 
description of the gas cooling process of the disc of M51a and 
adopt the local SF law in the model to calculate its star for¬ 
mation rate (SFR) using the cold gas surface density. Then, 
our model can predict the SFH, the chemical and colour 
evolution of M51a, with the help of stellar population syn¬ 
thesis (SPS) method. By comparing model predictions with 
the observed data, we can discuss the probable ranges for 
free parameters in the model and then know more about the 
main properties of the evolution and SFH of M51a. 

The paper is structured as follows. Section 2 describes 
the observed features of M51a disc, including the surface 
brightness, the cold gas content, the SFR and the metallic¬ 
ity. In Section 3 we present the main assumptions and in¬ 
gredients of our model in details. The comparisons between 
model predictions and observations and our main results are 
shown in Section 4. Our main conclusions are presented in 
Section 5. 


2 OBSERVATIONS 

M51a is one of th e closest (the distance: D ~ 8.0 Mpc; 
Walter et al.l 120081 ). face-on (inclination: i ~ 20°; iTullvl 
19741 1 SA(s)bc spiral galaxy a nd many of its b asic quantities 
have already been measured dHu et al.ll2013l l. A summary 


Table 1 . The main properties of M51a. 


Hubble type 


SA(s)bc 

i 

Distance 


8.0 Mpc 

2 

Inclination 


22° 

3 

M k 


-24.19 mag 

4 

Stellar mass 


~ 3.6 X 10 10 Mg 

5 

Hi mass 


~ (2.8 - 3.9) x 10 9 M 0 

5, 6 

H 2 mass 


~ (2.52 - 7.1) x 10 9 M 0 

5, 6 

/gas — -^Hgas/ (-^Hgas Af s tar) 

~ (0.13 - 0.23) 


SFR.puv+24/im 


3.125 M 0 yr —1 

5 

sfr Hq 


5.4 Mq yr —1 

7 

SFR R C-20cm 


2.56 Mq yr —1 

8 





Refs: (T) NED; (2) IWa 

ter et al. (200|); (3) Tully (1974) 

(4) 

Jarrett et al. (2003f); (5) 

Leroy 

et al, J2008T); (6) Miyamoto et al. 


d2014h : (71 iKennicutt et alj ll2003l ); (81 ISchuster et alj ll2007l l 


of the main properties of M51a is shown in Tabic [TJ most of 
which are adopted to constrain our model. 

In this Section, we summarize the current available ob¬ 
servations of M51a, especially the radial distributions along 
the disc, including the surface densities of gas mass and SFR, 
surface brightness and metallicity. 


2.1 Radial profiles of surface brightness 

The radial surface brightness profiles in multi-band of a spi¬ 
ral galaxy contain fossil information about the SFH along 
the galactic disc and provide strong constraints on the mod¬ 
els of galactic formation and evolution. 

The catalog data of the radial profiles of surface bright¬ 
ness for M51a in the ultraviolet bands (UV; Galaxy Evolu¬ 
tion Explorer, GALEX for short), optical bands (Sloan Dig¬ 
ital Sky Survey, SDSS for short), and near-infrared bands 
(NIR; Two-Micro n All-Sky Survey, 2M ASS for short) have 
been published bv lMunoz-Mateos et ali (120091 . MM09I here¬ 
after), but all these values were only corrected for fore¬ 
ground Galactic extinction and still in lack of internal ex¬ 
tinction c orrect ion. We adopt the method presented by 
iLin et ahl (120131) to estimate the internal extinction profile 
in the FUV band A fuv- Then, followin g the p rescriptions of 
iBoselli et al.l (120031) and ICortese et all ll2008h . we calculate 
the internal extinction in other wavelengths by assuming 
a given extinction law and a geometry for the distribution 
of stars and dust. The radial surface brightness profiles in 
FUV-, NUV-, u-, g-, r-, i-, z-, J-, H- and A-band are plot¬ 
ted in the left-hand side of Fig. [3] and Fig. |4j where the 
open circles show the observed profiles only corrected for 
the Galactic extinction and the filled ones also include a 
correction for the radial variation of internal extinction. 

The surface brightness in the A-band has another appli¬ 
cation. Since the 7\-band luminosity is most sensitive to the 
old stellar population and less effected by dust attenuation 
than other bands, we adopt t he A-band mas s -to-light ra- 
tio to be = 0.5 M 0 /Eq.k dBell et al.ll2003l : iLerov et ~akl 
l2008l) . and then utilize th e observed A'-band surface bright¬ 
ness profile (iMunoz-Mateos et al.l l2009h to estimate the 
present-day stellar mass surface density profile E*(r, t g ), 
where t g is the cosmic age and we set f g = 13.5 Gyr ac¬ 
cording to our adopted cosmology (see Section [3]). We will 
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use E*(r, t g ) derived here as our primary input parameter 
to constrain the total mass infalling to the disc. 

2.2 Radial profiles of cold gas and SFR 

During the past years, a number of data sets relating to 
the atomic and molecular gas surface density in M51a are 
becoming available. 

In spirals, the CO line emissions are usually used to esti¬ 
mate the molecular gas mass surface density Eh 2 • Imaging of 
molecular clouds has recently been carried out in M51a with 
the Berkeley-Illinois-Maryland Association (BIMA) i nterfer¬ 
omet er using the CO J = 1 — 0 line transition dHelfer et al.l 
l2003l l. Observations of the CO J = 2—1 emission in M51a 
were conducted with the Institute de RadioAstronomie Mil- 
limetrique (IRAM) 30m telescope using the 1 8 elem ent fo¬ 
cal plane heterodyne receiver array HERA dSchuster et ali 
120041') . Observations of CO J = 1 — 0 emission for M51a were 
also carried out using the 45-m telescope of the Nobeyama 
Radio Observatory (NRO) with th e 5 x 5—be am SIS hetero¬ 
dyne receiver array (BEARS'). ISchuster et al.l d2007l) derived 
Eh 2 with the help of the CO J = 2—1 line from HERA focal 
plane array, while others used the CO J = 1 — 0 line maps 
from BIMA and BEARS to obtain Sh „ (IKennicutt et aid 
120071 : iLerov et al.ll2008l : iMivamoto et alJl20l4 l. 

To derive the surface dens it y of n e utral atomic hydro - 
gen of M51a, Kennicutt et aid 120071 ') , ILerov et al.l d 20081 1 
and IMivamoto et al.l ~( 2014l l used Very Large Array (VLA) 
maps of the 21 cm line obtained as part of The H i 
Nearby Galax y Survey (THINGS) dWalter et al.l l2008tl . 
ISchuster et al.1 d2007T ) estimated the corresponding surface 
density from the large-scale distribution jrf the 21c m line 
of atomic hydrogen in M51a analyzed bv lRots et al.l dl990l l 
using the VLA. 

The sum of the surface density of atomic and molecu¬ 
lar gas is often called as the total gas surface density, i.e., 
Egas = 1-36(Eh 2 + Ehi), and the factor 1.36 is to include 
the contribution of helium. The observed radial distributions 
of molecular, atomic and total gas mass surface density are 
plotted in the right-hand si de of Fig . [5] and Fig. H] where the 
plott ed data are taken f rom Schuster et al.l (120071) (filled cir¬ 
cles), IKennicutt et al.l 1 20071') (filled triangles), ILerov et ah! 
J2008ll (asterisks) and IMivamoto et al.1 (120141 1 (filled dia¬ 
monds). 

With respect to the radial distribution of the SFR, sev¬ 
eral groups have already measured the SFR surface density, 
Esfr, in the disc of M51a using different tracers. The radio- 
continuum at 20 cm wavelength is used to est imate Esfr by 
ISchuster et al.l (120071 1. IKennicutt et al.l (l2007l l used Paa and 
a comb i nati on of 24 pm and Ha emission to obtain Esfr. 
ILerov et al . (2008) combined FU V and 24 pm maps to esti¬ 
mate Esfu . [Heesen et alj d2014l l obtained the radial distri¬ 
bution of the Esfr from A22 cm radio-continuum emission. 

The right fourth panel of Fig. [3] and Fig. [I] plots the 
radial profi les of SFR surface densi ty, Esfr- Th e data 
taken from ISchuster et al.l (120071 4 and iHeesen et al.l (120141 4 
are shown as filled circles (A20cm RC emission) and the 
filled pentagons (A22 cm RC emi ssion), respectively. The 
data obtained from ILerov et al.l d2008l l are represented 
by the the fille d asterisks (24 pm + FUV), while that 
from IKennicutt et al.l (120071 1 are denoted by filled triangles 
(24 pm + Ha). 


2.3 Radial profiles of metallicity 

The observed metallicity gradient is another important con¬ 
straint on the models of the galactic evolution. Except for 
hydrogen and heli um, oxygen is the most abundant element 
in the Universe (tKorotin et al.ll2014l f. and oxygen abundance 
is easily measured in Hu regions because of its bright emis¬ 
sion line. In practice, the oxygen abundance is sometimes 
used to represent the metallicity of the galaxy. 

The observed oxygen abundance gradients of Hu re¬ 
gions in M51a disc have been estimated by several authors. 
Based on analysis of a sample of 10 Hu regi ons w ith high- 
quality electron temperatures measurements, iBresolin et al.1 
(l2004ll derived an oxygen rad ial gradient of —0.021 ± 
0.011 dexkpc -1 . iMoustakas et akl (|2010|) adopted two cal¬ 
ibration methods to convert line strength to oxygen abun¬ 
dance and obtained two values of radial gradient. One is 
—0.038 ± 0.004 dexk p c" 1 by usin g the theoretical calibra¬ 
tion method of iKobulnickv fe Kewlevl (120041 . KK04 here¬ 
after) and the other is —0.024 ± 0.004 d exkpW 1 by ad opt¬ 


ing th e empirical calibr ation method o fl Pilvuein fe Thuaril 
l|2005l . PT05 hereafter). IPilvuein et al.l (12014 1 obtained an 
oxygen radial gradient of —0.0223 ± 0.0037 dexkpc -1 . 

The observed oxygen abundance profile of the M51a 
disc is shown in the right bottom panel of Fig. [3] and Fig. 
H The cycles re p resen t the observed data obtained from 


Moustakas et al.1 (hoiolb where the open and filled circles 


are derived by the PT05 and KK04 calibration method, re¬ 
specti vely. In the same panel, the data from IPilvuein et al.l 
(1201414 are shown as filled squares. 

However, the true situation of the oxygen abundance 
gradient in the M51a disc is still not fixed. The main un¬ 
certainty comes from the calibration methods used to derive 
chemical abundance. In addition, the data might suffer from 
lack of large enough samples. In any case, a larger and homo¬ 
geneous sample of Hu regions which spread the whole M51a 
disc is needed in order to have conclusive results about the 
real oxygen abundance gradients. 


3 THE MODEL 


Similar to previous models of the Milky Way 1 

Chang et al. 

19991; Chiappini et al.l 200ll: Matteucci et al. 

2003) and 

other nearby disc galaxies (Boissier & Prantzos 

2000; 

Molla & Dfaz 20051; Dalcanton 2007; Yin et al. 

2009; 

Munoz-Mateos et al.l 2011 

; Kang et al. 2012|; Chang et al. 

20l3: Robles-Valdez et al. 

20ll), we assume that the M51a 


disc is progressively built up by the infall of primordial gas 
(X = 0.7571, Y p = 0.2429, Z = 0) from its dark matter 
halo. The disc is basically assumed to be sheet-like and com¬ 
posed of a series of independently evolved rings with width 
500 pc, in the sense that no radial mass flows are considered. 

Throughout this paper we adopt the standard cold dark 
matter cosmology with Ho = 70kms _1 Mpc~ x , Dm = 0.3, 
and Ha = 0.7. Correspondingly, the cosmic time at z = 
0 is 13.5 Gyr. The details and essentials of our model are 
described as follows. 
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time (Gyr) 


Figure 1. The normalized gas infall rate (Upper panel) and the 
total infall gas surface density (Lower panel). Different line types 
correspond to various parameter groups: dashed lines (t p ,a) = 
(0.1 Gyr, 3.0 Gyr), dotted lines (t p ,cr) = (0.1 Gyr, 6.0 Gyr), dot- 
dashed lines (i p ,<r) = (15 Gyr, 6.0 Gyr), solid lines (£ p ,cr) = 
(15 Gyr, 3.0 Gyr). 


3.1 Gas infall rate 


We assume that there is a 1 Gyr time delay for the disc 
formation, which corresponds to z ~ 6 under the stan¬ 
dard cosmology. After that, the disc originates and grows 
by continuous primordial gas infall from the dark matter 
halo. We adopt a Gaussian f ormula of the gas infall rate 
from IChang et al.l (Il999l . 120121 ') . 

For the given radius r, the infall rate fi n (r,t) (in units 
of Mq pc -2 Gyr -1 ) is assumed to be: 




A( r ) -(t- to ) 2 / 2 ^ 
^2W 


(1) 


where t p is the infall-peak time and a is the full width at the 
half-maximum of the peak. The A(r) is a set of normalized 
quantities constrained by the stellar mass surface density at 
the present time £*(r, f g ). To explore how t p and a regulate 
the shape of the gas infall rate, Fig. |T] plots the normalized 
gas infall rate (upper panel) and the corresponding total in¬ 
fall gas surface density (lower panel), which is defined as how 
many cold gas have already been cooled to the given ring 
at time t and normalized to 1 at the present time. Different 
line types correspond to various parameter groups ( t p ,a ). 



Figure 2. The normalized quantities ,4(r). The dashed and solid 
lines are corresponding to two limiting cases of t p =0.1 Gyr and 
t p = 15 Gyr, respectively. 


We select t p = 0.1 Gyr and t p = 15 Gyr to represent the 
extreme case which corresponds to a time-decreasing infall 
rate and a time-increasing one, respectively. It can be seen 
from Fig.[Qthat the degeneracy of t p and a does exit in that 
the combination of small t p and large a results in a similar 
shape of mass accretion curve comparing that of large t p and 
small a. Since our main aim is to investigate the main trend 
of SFH for M51a but not to find the exact values of free pa¬ 
rameters, for the purpose of simplicity, we choose to reduce 
the number of free parameters in our model by treating t p 
as the free parameter and fixing a = 3 Gyr in the following 
section of this paper (IChang et al.|[201ol . 120121 ). 

We s h ould p oint out that, i n pre vious works by 
lHou_et all (l2000l ) . IChiappini et al.l ll200ll) and lYin et akl 
(l2009tl . an exponential form of gas infall rate with one free 
parameter is widely adopted in models of galactic chemical 
evolution. But, in these models, the gas infall rate always 
decreases with time and the most extreme case corresponds 
to a constant one. Therefore, we adopt a Gaussian formula 
of infall rate in this paper to include more possibilities, es¬ 
pecially the time-increasing gas infall rate. 

The normalized quantities A(r) in Equ. JT|) control how 
many cold gas finally cool to the disc as a function of ra¬ 
dius. In practice, we iteratively estimate A(r ) by requiring 
the model resulted stellar mass surf ace density at the present 
time is equal to i ts observed value (IChang et alJ feoiQ. 2012; 
[Kang et al.ll2012l) . In our calculation, after fixing the values 
of t p and cr, we give an initial value of A(r) and numeri¬ 
cally calculate the gas evolution and the increase of stellar 
mass. By comparing the resulted £«(r, t s ) with its observed 
value, we adjust the value of A(r) and repeat the calculation 
untill the resulted £„(r, t g ) is comparable to its observed 
value. Fig. [2] plots the resulted radial profile of A(r), where 
the dashed and solid lines are corresponding to two limiting 
cases of t p = 0.1 Gyr and t p = 15 Gyr, respectively. It shows 
that the A(r) does vary with t p . But we do not treat A(r) 
as free parameter in our model since it is roughly fixed for 
given t p . 

Regarding the stellar mass profile at the pr esent time, 
in our previous works such as IChang et al] d2012l) and 
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ll<ang et~al] (120121 1. we assumed that the stellar disc is a pure 
exponential disc and the value of E*(r, f g ) is constrained by 
the quantities of total stellar mass and radial scale-length, 
which are usually obtained from the observed radial profiles 
of surface brightness plusing the exponential disc assump¬ 
tion. In this paper, we regard E*(r, f g ) as a set of individual 
quantities and can be directly estimated from the observed 
radial profile of R'-band surface brightness by adopting a 
constant mass-to-light ratio (see Section 12.11 for details). 
Since this method includes much more spatial structures 
than the previous ones, it permits us to more reasonably 
calculate the multi-band radial profiles of surface brightness 
and compare with the observations in details. 


3.2 Star formation law 


SF process plays an important role in regulating the evolu¬ 
tion of galaxies. Unfortunately, the SF process is so compli¬ 
cated that its unde rlaying physical nature is still not com¬ 
pletely understood ( IKennicutt et al.ll2007l) . Most models of 
the galactic evolution still rely on the empirical SF law, 
which connects the local cold gas surface density to the SFR 
surface density, to describe how much cold gas converts into 


stellar mass in each t ime step dBoissieiLfc^ran tzos 
Cliiappini et alj : 200ll : lYin et aD 120091 ; Chang et al.l 
Kang et al.l 2012T U 


2000l; 


20121 : 


Based on the observed dat a of a samp l e of 9 7 nearby 
normal and star-burst galaxies, iKennicuttl (|199a l found a 
power-law correlation between the galaxy-averaged SFR sur¬ 
face density and the galaxy-averaged total gas surface den¬ 
sity, which was termed as the classical Schmidt-Kennicutt 
SF law and widely adopted in the studies of galactic forma¬ 
tion and evolution. But a more fundamental way is to use 
spatially resolved measurements of the SFR and gas sur¬ 
face density to examine the correlations between the ob¬ 
servables on a point-by-point basis within galaxies. Indeed, 
the studies of the spatially resolved SF law ha ve been ap¬ 
plied to several nearby spiral galaxies (IK ennicutt et al.ll2 007l : 
ILerov et al.ll2008l : iBigiel et~ak 120081 : ISchruba et al.ll20Lll) . In 
the case of M5 1 a, usi ng th e observed multiwavelen gth data, 
ISchuster et al.1 (l2007tl and IKennicutt et al.l (l2007li indepen¬ 
dently investi gated the spat ially resolved SF law along the 
disc of M51a. IKennicutt et al] (J2007]) found that the resolved 
SFR versus gas surface density relation is well represented 
by a Schmidt power law, which is similar in form and dis¬ 
persion to the disc-average Schmidt-Kennicutt law. 

In our model, we ado pt the SF law of M51a presented 
by IKennicutt et ahl (1200711 : 

'h(r, t) = 0.048E gas (r, t) 1 ' 56 , (2) 

where >3(in units of Mq pc -2 Gyr -1 ) is the SFR sur¬ 
face density, and E gas (r,t) (in units of Mq pc -2 ) is the 
gas surface density, which is the sum of the surface den¬ 
sity of the atomic and molecular gas. We assume that 
E gas (r,t) = 1.36(Eh 2 (r,t) + S H i(r, t)), where E H2 (r,t)and 
Emi(r, t ) are the surface density of atomic and molecular hy¬ 
drogen, and the factor 1.36 is to reflect the contribution of 
helium. 

Regarding the ratio of molecular-to-atomic gas sur- 
face density A m oi(r, £), we ado pt the formu l a pre sented by 
[Blitz fc Rosolowskvt (|2006i l and ILerov et al.l ((2008h , 

A mo i(r,f) = Eh 2 (r, f)/E Hl (r,t) = [P h (r,t)/P 0 ] ap , (3) 


where Ph(r,t) is the mid-plane pressure of the interstellar 
medium (ISM), Po and ap are constants derived from the 
observations. We adopt Pp/k = 1.92 x 1 0 4 cm -3 K and ap = 
0.87 derived bv IHitschfeld et ~akl (l2009l l for M51a. 

The mid-pl ane pressure of the ISM in disc gala xies can 
be expressed as (lElmegreenlll989l : ILerov et alJl2008h : 


Ph(v,i) — — GEgas (r, t) 


E gas (r,t) + ^E„ (r, t) 
c* 


(4) 


where G is the gravitational constant, and c gas and c* are the 
(vertical) velocity dispersions of gas and stars, respectively. 
Observations reveal that c gas is a constant along the disc and 
we adopt e gas = llkms -1 dOstriker et al.ll2010ll . but c» is 
estimated as c* = y/ 7rGzoE* (r), where 20 is th e scale-heig ht 
of th e disc and we adopt 20 = 1.0 kpc (IHitschfeld et al.l 
l2009ll . 


3.3 Other ingredients and basic equations 

We also take into account the contribution of the gas out¬ 
flow proc ess i n our model. According to the mass-dependent 
model of IChang et alJ (hold ). we assume that the gas out¬ 
flow rate is proportional to the SFR and the coefficient is 
set to be & 0 ut = 0.004. It should be emphasized that our 
final results is not sensitive to the variation of 6 ou t since the 
outflow process plays a relatively small role in the chemical 
evolution of such a massive galaxy as M51a (the stellar mass 
of M51a is about M* = 10 10 ’ 7 Mq). 

W e adopt the updated SPS model of lBruzual fe Charlotl 
d2003ll (i.e., CB07) wit h the stel l ar ini tial mass function 
(IMF) being taken from IChabrierl d2003h in our work. The 
lower and upper mass limits are adopted to be 0.1 Mq and 
100 Mq, respectively. 

Regarding the chemical evolution of the disc of M51a, 
both the instantaneous-recycling approximation and the in¬ 
stantaneous mixing of the gas with ejecta are assumed, that 
is, the gas in a fixed ring is characterized by a unique com¬ 
position at each epoch of time. We take the classica l set of 
equations of galactic chemical evolution from lTinslevI (1 1980h : 

= fin(r,t) - /out(r,t), (5) 


d[Etot 


rfpWr.t)] = -(l-A)*(r,f) + /inM)-/outM), (6) 


d[Z(r,f)S gaB (r,t)] = y(1 _ m{r , t) _ Z { r ,t)G - R)*(r,t) 

+Z in fin(r,t) - Z ou t(r, t)f ou t(r, t). (7) 

where E to t (r, t) is the total (star + gas) mass surface density. 
Z(r,t) is the metallicity in the ring centered at galactocen- 
tric distance r at evolution time t. R is the return fraction 
and we set R = 0.3 according to the adopted IMF. y is 
the stellar yield and we set y = 1 Zq (IChang et al.l l20ld : 
iKang et al.ll2012l ') . Zi n is the metallicity of the infallmg gas 
and we assume the infalling gas is primordial, that is Z i n = 0. 
Z ou t(r,t) is the metallicity of the outflowing gas and we as¬ 
sume that the outflow gas has the same metallicity as that of 
ISM,e.g., Z out (r,t) = Z(r,t) (IChang et al.l2Qlol : lKang et al.l 

l2012h . 

We emphasize that, under the condition that we have 
already fixed er and A(r), there is only one free parameter 
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Figure 3. Comparisons of the radial profiles between model predictions and observations. The dashed and dotted lines are corre¬ 
sponding to the model results adopting t p = 0.1 Gyr and t p = 15Gyr, respectively. In the left-hand panels, the open circles show 
the observed profiles only corrected for the Galactic extinction and the filled ones also include a correction for the radial variation of 
internal extinction. The observed dat a of radial distr i butio ns of molecular, atomic, and total gas mass surface dens ity shown in panels 
of the right-ha nd side are taken from Schuster et alj (120071 ) (filled circles). iKennicutt et all ([2007 ) (fill ed triangles ). I Leroy et al.l (120081) 
( aste risks) and [ Miyam oto et al] (12014 ) (fille d diamonds) , and the SFR surface density obtained fromfechuster et al.l 120071) (filled cir¬ 
cles), [Heese^et^ljli^SMfiiiGd pentagons). I Lerov et al.1 (|2008l ) faste risks) and [K ennic utt et al.1 (|2007l ) (filled triangles) are displayed in 
the right-f ourth panel. The observ ed oxygen abundances taken from [Pilvugin et al.l (|2014|) are shown as filled squares, while the data 
taken from lMoustakas et al.1 (|20ld ) are plotted as open and filled circles corresponding to be derived by the PT05 and KK04 calibration 
method, respectively. 
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Table 2. The model predictions of the total quantities of M51a 
with different parameters 


tp 

(Gyr) 

M k 

(mag) 

Mh 2 

GO 9 M 0 ) 

m Hi 

(10 9 M 0 ) 

is 3 

0 Pd 

i-i 

1 

/gas 

0.1 

-24.17 

1.86 

2.16 

1.04 

0.07 

7.0 

-24.84 

6.92 

2.96 

5.02 

0.17 

15.0 

-26.08 

38.22 

3.21 

49.75 

0.46 


left in our model, the infall-peak time t p , which regulates the 
shape of gas accretion history and then largely influences the 
main properties of SFH along the disc and then the evolu¬ 
tion of M51a. Generally speaking, t p —> 0 is corresponding 
to a time-declining infall rate, while t p —» oo is correspond¬ 
ing to a time-increasing gas infall rate. By comparing model 
predictions with the observed data, we can discuss the prob¬ 
able range of the free parameter in the model and then know 
more about the main properties of the evolution and SFH 
of M51a. 


4 MODEL RESULTS VERSUS OBSERVATIONS 

In this section, we present our results step by step. Firstly, 
we investigate the influence of the free parameter on our 
model predictions and present the viable model of the evo¬ 
lution of M51a. Furthermore, we explore the properties of 
stellar population along the disc and compare the growth 
history of M51a with that of UGC 8802 and the Milky Way. 


4.1 Radial profiles 

To explore the influence of free parameter on model results, 
we firstly consider two limiting cases of t p = 0.1 Gyr and 
t p = 15 Gyr and present the comparison between model pre¬ 
dictions of the radial profiles and the observations in Fig. [3j 
The dashed and dotted lines in Fig.[3]denote the predictions 
of the models adopting t p = 0.1 Gyr and t p = 15 Gyr, re¬ 
spectively. The left-hand side of Fig. [3] displays the radial 
multi-bands surface brightness profiles from FUV-band to 
R'-band, where the open circles show the observed profiles 
only corrected for the Galactic extinction and the filled ones 
also include a correction for the radial variation of internal 
extinction. The right-hand side of Fig. [3] shows the H 2 , HI 
and total gas mass surface density, SFR surface density, and 
oxygen abundance radial profiles. The details of the observed 
data are presented in Section [2] 

Large differences between dashed and dotted lines in 
Fig. [3] shows that the model predictions are very sensitive 
to the adopted t p . The case of t p — 0.1 Gyr (dashed lines) 
corresponds to a time-decreasing gas-infall rate that most of 
the gas has been accreted to the disc in the early period of its 
history, while that of t p = 15 Gyr (dotted lines) corresponds 
to a time-increasing gas infall-rate that a large fraction of 
cold gas is still infalling to the disc at the present time. It 
can be seen that the model adopting the earlier infall-peak 
time (dashed lines) predicts less luminous surface bright¬ 
ness, lower molecular gas surface density, lower SFR sur¬ 
face density and higher gas-phase oxygen abundance than 


that adopting a later infall-peak time (dotted lines). This 
is mainly due to the fact that, in our model, the setting of 
earlier infall-peak time corresponds to a faster gas accretion 
and higher SF process in the earlier stage of its history and 
then leads to the older stellar population age, higher gas- 
phase metallicity and lower cold gas content at the present 
day. 


Interestingly, we can see from Fig. [3] that the area be¬ 
tween the dashed and dotted lines brackets almost the whole 
region of the observations, which indicates that it is possible 
to construct a model that can reproduce the main features of 
the observations of the M51a disc. In fact, these two extreme 
cases (t p = 0.1 Gyr and t p = 15 Gyr) are bracketing the 
possible range of our model results, from where we can see 
the variation trends of the model predictions with different 
value of t p . After a series of calculations and comparisons, 
we obtain a viable model by adopting the infall-peak time 
tp/Gyr = 7.0 and plot its results as solid lines in Fig|4] The 
notation of observed data is the same as that of Fig[3] We 
can see from Fig. U that the solid lines can basically repro¬ 
duce most of the observational data, which indicates that 
our viable model includes and describes reasonably the key 
ingredients of the main processes that regulate the formation 
and evolution of M51a. 

To further understand the effect of t p on the evolu¬ 
tion of M51a, Table [2] summarizes the model predictions 
of the present total quantities of M51a, including the abso¬ 
lute K-band magnitude Mk, the molecular gas mass Mh 2 , 
the atomic gas mass Mm, SFR and the gas fraction / gaa . 
Three different values of t p (0.1, 7.0,15Gyr) are considered, 
which corresponds to the dashed, dotted lines in Fig. [3] and 
solid lines in Fig(4] It can be seen from Table [2] that the 
model adopting a later infall-peak time predicts higher SFR 
and higher gas fraction. The predictions of the viable model 
adopting f p /Gyr = 7.0 are in fairy agreement with the cor¬ 
responding physical properties in Table |T] considering the 
observed uncertainties. 

However, it can be seen from Fig[I] that, contrary to 
the smoothness of the model predictions, the observed ra¬ 
dial profiles of cold gas and SFR surface density show peaks 
at Rg ~ 6 kpc and the location of these peaks coincide with 
the location of the most a ctive SF region in the disc of M51a 
llScheepmaker et al.ll200sh . In addition, bumps obviously ap¬ 
pear on the surface brightness radial profiles in the FUV-, 
NUV- and u-band, and the intensity of the bump decreases 
gradually with the increase of the wavelength, which sug¬ 
gests that active SF process may take place in recent years 
around there. 

Indeed, both kinematic and hydrodynamic modeling 
suggest that single or multiple encounte rs between these two 
galax i es are occurred 300-500 Myr ago (Salo &; Laurikainenl 
l2000l : iTheis fe Spinnekeil 120031 : iDobbs et alJl2Qloh . Studies 
of stellar populations of M51a also revealed that a b urst of 
SF occurred 340—500 Myr ago llMentuch Coope r et al .11 201 i ; 

I Tikhonov et al.l 120091') . Moreover. iMivamoto et al.l (12014 ) 
found that the local maximum at Rg ~ 6 kpc on the ra¬ 
dial profiles of cold gas in M51a could be caused by the 
kinks or fractures of the spiral arms by the interaction with 
its companion galaxy. 

Correspondingly, we present a ’toy’ model to imitate the 
influence of the interaction between M51a and its compan¬ 
ion by allowing an additional cold gas infall (with constant 
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Figure 4. Comparisons of the model predictions with the observations. The solid lines plot the viable model results, while the dashed 
lines represent the ’toy’ model results. The notation of the observational data is the same as that of Fig. [3] 


rate 8.31 Mgyr -1 ) to the disc occurred recently (during the 
period 13.0 ^ f/Gyr ^ 13.2) round the region Rq ~ 6 kpc. 
Other ingredients of the model are the same as those of the 
viable model. The results of the ’toy’ model are shown by 
dashed lines in Fig. U It can be seen that the freshly in¬ 
falling cold gas increases the mass surface density of both 
the molecular gas and the SFR, while has little effect on the 
atomic surface density and the gas-phase oxygen abundance. 
Moreover, the additional gas infall results in more luminous 
surface brightness in FUV-, NUV- and u-band, while the sur¬ 
face brightness in optical and near infrared bands are less 


effected. Our results suggest that the molecular gas surface 
density, the SFR and the UV-band luminosity are impor¬ 
tant quantities to explore the effects of recent interaction on 
galactic SF process. 

4.2 Stellar populations along the disc 

To further explore the properties of stellar population along 
the disc of M51a, we select four regions out and plot their 
viable model predictions of the SFHs and growth curves 
of stellar masses in Fig. [5] including 1-2kpc (dashed line), 















The Evolution of Interacting Spiral Galaxy NGC519f 9 



Figure 5. Evolution history of the M51a disc predicted by the 
viable model. Upper panel: the model predicted time evolution 
of SFR for four main spatial components (1-2 kpc, 4-5 kpc, 7- 
8 kpc, 10-11 kpc) for M51a disc. Lower panel: relative stellar mass 
growth of the corresponding four spatial components as the upper 
panel. Both the SFRs and stellar masses are normalized by their 
present-day values. The horizontal dash-dotted line in the lower 
panel remarks when each component achieves 50% of its final 
value. 


4-5kpc (dotted line), 7-8kpc (dash-dotted line), 10-11 kpc 
(solid line). Correspondingly, Fig. [6] displays their metal- 
licity distribution function (MDF) predicted by the viable 
model. Note that both the SFRs and stellar masses in Fig. 
[5] are normalized by their present-day values, and the hori¬ 
zontal dash-dotted line in the lower panel of Fig. [5] denotes 
the position that each component achieves 50% of its final 
value. 

Fig. [5] shows that the SFRs in the whole disc are very 
low at early stage and increase gradually due to the increase 
of gas infall rate. When the consumption rate of cold gas to 
form new stars is roughly balanced by the gas infall rate, the 
SFR reaches its peak and then drops down to its present- 
day value. It is interesting to see that the peak of the SFH 
moves gradually to later time from the inner to outer parts 
of disc. Moreover, Fig. [G] reveals that the inner disc contains 



Figure 6. The viable model predicted metallicity distribution 
function for four main spatial regions (1-2 kpc, 4-5 kpc, 7-8 kpc, 
and 10-11 kpc) of M51a disc. 

a higher percentage of metal-rich stars than the outer disc. 
In other words, our model predictions show evidences that 
the disc forms inside-out. We will discuss this point again in 
next subsection. 


4.3 Comparisons with the Milky Way and 
UGC 8802 

The Milky Way and UGC 8802 are two disc galaxies con¬ 
taining comparable stellar mass with M51a, and we have 
already investigated SFHs of th e two galaxies in previous 
works dChang et al.l 1 19991 . 120121 1 . In this section, we will 
compare the growth history of M51a predicted by the vi¬ 
able model with that of the M ilky W ay and UGC 8802. We 
referred to IChang et al.l 1 19991 . 12OI2I) for a more indepth de¬ 
scription of the models. Here we only repeat main ingredi¬ 
ents of the models such as the gas infall rate and the SF law, 
and present a summary of the models in Table [3] 

In our previous models of the Milky Way and 
UGC 8802, an Gaussian formula of gas infall rate has been 
adopted. The input condition of UGC 8802 to constrain how 
much gas infalling to the disc is the stellar mass surface 
density at the present day, which is assumed to be pure ex¬ 
ponential S,(r ,t g ) = 93.7e~ r / rd - u , with radial scale-length 
rd,u = 5.8 kpc IlChang et al.l[201^ ). In the case of the Milky 
Way, the input parameter to constrain the gas infall rate 
is the total (gas+star) mass surface density at the present 
day, i.e., Stot(r,t g ) = 55.0e _ ^ _7 '®^ rd ' mw , where the so¬ 
lar Galactic radius is set to b e r^ = 8.5 k pc a nd the disc 
scale-length is ra.mw = 2.7kpc dChang et alJll99Sh . Regard¬ 
ing the SF law, the model of UGC 8802 adopted the one 
that the SFR is proportional to the molecular hydrogen sur¬ 
face density fk(r, t ) = E mo i(r, t)/ 0.77, while the model of the 
Milky Way adopted a radial-dependent Schmidt SF law, i.e., 
’F(r, t) = 0.85 Eg a g(r,t)(r/rd, mw ) _1 , to be consistent with 
the observed steep oxygen gradient along the disc. 

In Fig. [T] we compared model predictions of the growth 
history of M51a (solid lines), UGC8802 (dashed lines) and 
the Milky Way (dotted lines). The upper panels plot the evo¬ 
lution of stellar mass M, (left panel) and atomic hydrogen 
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Table 3. A summary of the viable models for M51a, UGC8802 and the Milky Way. 


Individual 

M51a 

UGC 8802 

Milky Way 

Input properties 

MK (0) 

= O.5M 0 /L 0 , k 

E*(r,tg) = 93.7 e -r '/ rd ’ u M 
r d,u = 5.8 kpc 

StotMg) = 55.0e - ^’ r 0)/ r a,mw (d) 
r 0 = 8.5 kpc and r < j jInw = 2.7kpc 

SF law 

V(r,t) = 0.048 EgasM) 1 ' 66 

^(r,t) = £ mo i(r,t)/0.77 

y(r,t) = 0.315 Egasiy, i) 1 ' 4 (yAd, mw ) — 1 

Infal-peak time t p (Gyr) 

7.0 

1.5r/rd lU + 5.0 

2.0r/r < j jMW + 2.0 

Notes. 


The iC-band mass-to-light ratio taken fro mlBe ll et al.| (2005) and lLerov et alj <(21 


rd )U is the disc scale-length of UGC 8802 \Chang et al.l 


20121 ) . 


r© is the solar Galactic radius, and mw is the disc scale-length of the Milky Way from [Chang et all <119991) . 




time (Gyr) 


0 3 6 9 12 

time (Gyr) 


Figure 7. The growth history of M51a, UGC 8802 and the Milky Way predicted by their individual viable models. The upper panels plot 
the model predicted time evolution of M* (left panel) and Mh x (right panel) of the three galaxies, The growth history of Mh 2 and the 
SFHS are displayed in the middle-left and middle-right panels, respectively. The model predicted evolution of the half mass size R 50 ,* 
is displayed in the bottom-left panel. Each quantity is normalized to its value at the present day. The different line types correspond to 
different galaxies: solid lines for M51a, dashed lines for UGC 8802, and dotted lines for the Milky Way. 
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mass Mh t (right panel). The growth history of molecular 
hydrogen mass Mh 2 an d the SFHs of these three galaxies 
are displayed in the middle-left and middle-right panels, re¬ 
spectively. The bottom left panel shows the evolution of the 
half mass size -R 50 ,*, which is defined as the radius at which 
half of the total stellar mass is contained. Each quantity in 
Fig. m is normalized to its value at the present day. 

From Fig. [7] we can see that, at early stage of evolution, 
most of the infalling cold gas is in the form of atomic gas. 
After atomic gas accumulates significantly enough, molec¬ 
ular gas begins to form and then the SFR speeds up. The 
middle panels of Fig. 0 show that the shape of the SFH is 
very similar to the evolution curve of the molecular gas, no 
matter whether the adopted SF law assumes the SFR pro¬ 
portional to the molecular gas surface density or not. This 
point suggests that the density of molecular gas may be the 
main driver of SF process, but this question is far from set¬ 
tled. 

The growth curves of the mass of atomic gas show that, 
both M51a and the Milky Way reach their peaks around 
4 ~ 7 Gyr ago and then slowly drop down to their present- 
day values, while the atomic mass of UGC 8802 increases al¬ 
most steadily and seems to just approach to its peak. Indeed, 
the outstanding point of UGC 8802 is its extremely high neu¬ 
tral gas content compared to other disc gala xies in this stel ¬ 
lar m ass range (jGarcia-Anpadoo eUah| [2009: Cati nella et al.1 
l20ldf . According to the model of lChang et al . I d2012h . a late 
infall-peak time is necessary to explain the high neutral gas 
content of UGC 8802, which implies that UGC 8802 is a 
young galaxy and still active in SF processes. This point 
can also be seen from both the upper-left panel and the 
middle-right panel of Fig. 0 It is shown that half of the to¬ 
tal stellar mass of M51a has been assembled during the last 
~5.0 Gyr, while that for the Milky Way and UGC 8802 is 
around ~6.5 Gyr and ~4.0 Gyr ago, respectively. I 11 other 
words, our results suggest that the stellar population of 
M51a is younger than that of the Milky Way, but older than 
that of the gas-rich disc galaxy UGC 8802. 

We also compare the evolution of the half-mass size 
R 50 ,* among these three galaxies in the bottom left panel 
of Fig. [7] Since the models of the Milky Way and UGC 
8802 both assume that the infall-peak time increases with 
radius (see Tabled, it is natural to see that both the dotted 
and dashed lines increase with time. In the case of M51a, 
although the model adopt a constant infall-peak time, the 
model predictions also support the disc inside-out forma¬ 
tion scenario. This trend also reinforces the recent analysis 
of the SFH of CALIFA galaxies, which shows that galax¬ 
ies more massive than IQ 10 Mp grow inside-out llPerez et al.l 
120131 ; IGonzalez Delgado et all 120141 '). In fact, this inside- 
out formation mechanism has already successfully ap¬ 
plied to previous models of formation and evolution of 


disc galaxies ( Chang et al. 1999, 2012 

; Boissier & Prantzod 

20001; Hou et al. 200(1 ChiaDDini et al. 

2 OOII; Yin et al]|2009j; 

Kane et al. 2012f) and supported by the observations for 


nearb y disc galaxies jMunoz-Mateos et~ al. 2007; Perez et al.l 
120131') . 


5 SUMMARY 

M51a is a grand-design spiral galaxy and well known to be 
interacting with its companion. In this paper, we focus on 
investigating the SFHs of M51a by constructing a param¬ 
eterized model and comparing model predictions with the 
observed data, especially the radial distributions of cold gas 
surface density, metallicity, and the radial profiles of surface 
brightness in multi-bands. Our main results can be summa¬ 
rized as follows: 

(i) Our results show that the model predictions are very 
sensitive to the adopted infall-peak time t p . A late infall- 
peak time t p results in more luminous surface brightness, 
low metallicity, high gas and SFR surface densities. 

(ii) The model adopting a constant infall-peak time t p = 
7.0Gyr can nicely reproduce most of the observed con¬ 
straints of M51a. We also introduce a ’toy’ model, which 
allows a small amount of additional cold gas infall to the 
disc in recent time, to imitate the influence of the inter¬ 
action between M51a and its companion. Our results show 
that the additional gas infall could be attributed to the ob¬ 
served small bumps on the radial profiles of cold gas surface 
density, SFR surface density and surface brightness in the 
UV-band. 

(iii) Although we adopt a constant infall-peak time, our 
model predictions still show evidence that th e disc f orms 
inside -out. This is consistent with the result s of lPerez et alJ 
(l2013h and IGonzalez Delgado et al] (120141 ') that galaxies 
more massive than 10 10 Mq grow inside-out. 

(iv) We also compared the model predicted growth histo¬ 
ries of M51a, the Milky Way and UGC 8802. We find that 
the mean stellar age of M51a is younger than that of the 
Milky Way, but older than that of the gas-rich disc galaxy 
UGC 8802. Our results also show that half of the total stel¬ 
lar mass of M51a may has been assembled during the last 
~5.0 Gyr and SF may still proceed actively in the disc of 
M51a. 

We emphasize that here we only present a parameterized 
model. Through the comparison between model predictions 
and observations, our ultimate goal is to present our story of 
the main properties of the formation and evolution of M51a. 
Although the accurate value of t p in the viable model is 
not unique, the main conclusions of our results are robust. 
However, since M51a belongs to a dynamically complicted 
system and maintains grand-designed spiral structure, it is a 
long way to go to clearly understand its origin and evolution. 
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